Random-Manifold to Random-Periodic Depinning of an Elastic Interface 
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We study numerically the depinning transition of driven elastic interfaces in a random-periodic 
medium with localized periodic-correlation peaks in the direction of motion. The analysis of the 
moving interface geometry reveals the existence of several characteristic lengths separating different 
length-scale regimes of roughness. We determine the scaling behavior of these lengths as a function of 
the velocity, temperature, driving force, and transverse periodicity. A dynamical roughness diagram 
is thus obtained which contains, at small length scales, the critical and fast-flow regimes typical of 
the random-manifold (or domain wall) depinning, and at large length-scales, the critical and fast-flow 
regimes typical of the random-periodic (or charge-density wave) depinning. From the study of the 
equilibrium geometry we are also able to infer the roughness diagram in the creep regime, extending 
the depinning roughness diagram below threshold. Our results are relevant for understanding the 
geometry at depinning of arrays of elastically coupled thin manifolds in a disordered medium such 
as driven particle chains or vortex-line planar arrays. They also allow to properly control the effect 
of transverse periodic boundary conditions in large-scale simulations of driven disordered interfaces. 

PACS numbers: 75.60.Ch, 64.60. Ht, 05.70.Ln 



I. INTRODUCTION 



The dynamics of clastic manifolds in disordered me- 
dia have been widely studied in relation with the 
physical properties of many systems. Magnctioir— or 
ferroelectrio^^— domain walls, contact lines of liquid 
menisci^i^ fluid invasion fronts in porous mediaj^iii and 
fractureai^rJ^ can be modelled as elastic interfaces or 
lines. Such lines would be flat if they were not under 
the usually unavoidable action of quenched disorder. Pe- 
riodic systems like charge density waves^^ vortex lattices 



in type II superconductors 
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or Wigner crystalsai can 



be also modelled as elastic manifolds embedded in ran- 
dom environments. This manifold is described by the dis- 
placements around the perfect periodic lattice that would 
exist in the absence of disorder. No matter how weak 
is the disorder^SS in all these systems the competition 
between elasticity and disorder gives rise to rough struc- 
tures and complex collective pinning phenomena with in- 
teresting universal features. 

Of special interest is the response of this kind of sys- 
tems to an external uniform field, able to drive the elas- 
tic manifold in a given direction. Concrete examples are 
applied magnetic fields on magnetic domain walls, ap- 
plied electrical fields on ferroelectric domain walls, fluid 
pressure on contact lines, tension on fractures, electrical 
currents on vortex lattices in superconductors, and elec- 
trical fields on charge density waves and Wigner crystals. 
Indeed; such a probe would be rather trivial if not be- 
cause of the presence of quenched impurities: Disorder 
breaks the translation symmetry (though not in a statis- 
tical sense), making the otherwise uniform displacement 
of the manifold a complicated process involving many 
degrees of freedom. Whether the elastic bonds of the 
manifold break or support the tearing produced by the 
disorder, the resulting fiow can be plastic or elastic re- 



spectively, and in both cases a rich disorder-induced out 
of equilibrium phenomena can emerge. To understand 
these phenomena it was shown to be more convenient to 
start by restricting the study to the more tractable elastic 
flow case. 

Elastic depinning is one of the most prominent and bet- 
ter understood examples of collective pinning dynamic 
phenomenai^i^i At zero temperature the external field 
must overcome a finite threshold fc in order to force the 
pinned system to acquire a finite steady-state velocity 
V. Below the depinning threshold a finite velocity is 
only possible at a finite temperature by thermal acti- 
vation, due to the presence of many metastable states 
separated by energetic barriers. These barriers tend to 
diverge when decreasing the drive in the so-called creep 
regime )^^'^^ strongly impeding the motion at low driving 
forces, and tend to vanish at fc giving place to a ther- 
mal rounding^ of the depinning transition. These col- 
lective transport phenomena are experimentally relevant 
since a finite velocity in this kind of systems correspond 
to physical quantities (magnetization, or polarization for 
domain walls, voltage for superconductors, current for 
CDW) that can be readily measured. 

From the statistical physics point of view the most re- 
markable feature of the far from equilibrium steady-state 
motion near the depinning threshold fc at zero tempera- 
ture is the existence of a well defined non-trivial critical 
behavior. Just above the threshold the motion is jerky, 
characterized by forwardly moving avalanches of a typical 
size ^ ^ {f — fcY and width w ^ produced at a typ- 
ical rate r ^ yielding a mean velocity v ^ [f — fc)^, 
with /3 — j/((^ ~ z). v, z, ^ are non-trivial characteristic 
exponents. These observations led to the fruitful anal- 
ogy of the depinning transition with standard equilibrium 
critical phenomena, with v playing the role of the or- 
der parameter and / the role of the control parameter 1^ 
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This analogy motivated an outburst of analytical and nu- 
merical work devoted to determine the value of critical 
exponents for different universality classes ;22i22ri^ and 
to develop powerful analytical^"— and numerica l^"'^^" — 
methods to obtain them. From the numerical viewpoint 
such a study requires a precise determination of the crit- 
ical threshold 

For standard equilibrium phase transitions the low- 
temperature phase can be characterized by equilibrium 
correlation lengths separating the critical-looking short 
length-scales from the low-temperature fixed-point dom- 
inated large length-scales. For the depinning transition 
it was shown that ^ also admits an analogous purely ge- 
ometric interpretation as a crossover length in the av- 
erage steady-state roughness of the ('ordered') moving 
w > phase. The length ^ separates the regime of crit- 
ical roughness at short length scales (i.e. with a rough- 
ness exponent of the critical configuration at /c) from 
the fast-flow roughness observed at large length-scales 
(i.e. with a roughness exponent identical to the strongly 
driven interface, / ^ /c)'^'^'^ The steady-state geom- 
etry thus contains information of the velocity and there 
is no need to observe the transient correlated process of 
an avalanche. More recently, however, the analysis of 
the low-temperature averaged steady-state geometry has 
shown that no divergent steady-state correlation length- 
scale exists approaching the threshold from below, thus 
breaking the naive analogy with standard phase tran- 
sitions, where two divergent length-scales are expected 
above and below the critical point 4^1^ 

Elastic depinning universality classes were shown to 
depend on the dimension of the embedding space D, the 
dimension d of the manifold or the number N = D — d 
of displacement components of the manifold, the na- 
ture of the elastic interactions r^iiSi^iS the anisotropy 
of the mediun*^ and the nature of microscopic disorder 
correlations i^^i^ Considering for simplicity the case of d- 
dimensional directed manifolds with A*" = 1 living in an 
isotropic uncorrelated disordered medium it is convenient 
to distinguish between two prominent groups, according 
to the correlations of the effective pinning force Fp{iL,Y). 
This pinning force acts on the manifold displacement field 
u{y), which measures the distance between the distorted 
and the perfectly flat manifold at the labelling point r. 
On one hand the pinning force on interfaces such as do- 
main walls or contact lines in random potentials usually 
display short-range correlations reflecting the fact that 
the interface sees a completely different disorder after 
shifting it a distance bigger than a certain characteristic 
finite width rf = max[i(j, tq], where w is the domain wall- 
width and ro the assumed finite correlation length of the 
disorder potential. We use Random-Manifold (RM) to 
denote this group, and we do not make distinction be- 
tween the Random-Bond and Random-Field type of dis- 
order since at depinning, unlike statics, they are known 
to merge into a single classi^i^ Interfaces in periodic 
potentials or periodic condensates such as charge den- 
sity waves or periodic chains of elastically coupled ob- 



jects on the other hand display an effective pinning force 
with periodic correlations with a period M representing 
the period of microscopic potential in the first case and 
the lattice spacing in the second case. We use Random- 
Periodic (RP) to denote this group. 

For short-range correlated isotropic disorder, the N ~ 
1 RM and RP classes have been traditionally studied, 
both numerically and analytically, using two paradig- 
matic models of disorder. While for modelling the large- 
scale dynamic behavior of a non-periodic system it is 
enough to use any uncorrelated potential with range r^, 
for modelling the periodic system the random-phase co- 
sine potential have been traditionally chosen, thus forc- 
ing M ~ rf. Although this is a good approxima- 
tion for charge density waves (cf. Fukuyama-Lee-Rice 
model^i^) , this kind of modelling does not permit how- 
ever to study the interesting situation that can appear 
in different periodic systems for which the periodicity is 
much larger than the short-range correlation length of the 
disorder correlator, i.e. r/ ^ M. Indeed, when the auto- 
correlation of the pinning force is periodic and displays 
sharply localized peaks this physical situation, mostly an- 
alyzed for two component {N = 2) displacement fields, 
was shown to be relevant for describing the statics of 
Wigner crystals^ii^ or vortex lattices j^i^ where the lat- 
tice spacing qq can be made much larger than the vor- 
tex core size or coherence length ^ (cyclotron radius for 
Wigner crystals) by simply tuning an external magnetic 
field. In these cases the length-scale separation is re- 
sponsible for the so-called RM regime of roughness. This 
regime occurs at intermediate length-scales, before the 
system asymptotically reaches the so-called Bragg-glass 
or RP regime. Because the intermediate RM regime can 
span a wide range of lengths r^i^ it can affect the stat- 
ics and dynamics properties of this kind of systems and 
thus can be experimentally observed. We can therefore 
expect additional geometrical crossovers around the de- 
pinning transition in these systems. From a numerical 
point of view, the effect of a periodicity M ^ rf has 
been already analyzed in the critical depinning force dis- 
tribution in Ref . l63l . 

Here we present a study of the finite velocity dynam- 
ics of a simple RP system which includes localized pe- 
riodic correlation peaks with controlled periodicity M, 
yielding an interesting multiscale behavior around de- 
pinning. Our main result is a geometrical dynamical 
roughness diagram which contains, at small length scales, 
the critical and fast-flow regimes typical of the RM (or 
"magnetic domain wall") depinning, and at large length- 
scales, the critical and fast-flow regimes typical of the 
RP (or "charge-density wave" ) depinning. We argue that 
our results arc qualitatively valid for the family of one- 
component periodic systems with localized correlations 
peaks, such as chains of elastically coupled thin inter- 
faces. We compare in particular a driven chain of inter- 
acting particles in a one-dimensional disordered potential 
with an elastic line in disordered potential with periodic 
correlations at a larger scale. Our results are particu- 
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larly relevant for properly controlling and interpreting 
the effect of periodic boundaries conditions in large-scale 
driven interface simulations. 



Outline of the paper 

The paper is organized as follows. In Sec. [TT] we de- 
scribe the general class of random periodic systems with 
localized periodic correlation peaks for which we argue 
our general results apply. Then, Sec. UTT] presents the 
general properties of the structure factor, which will be 
used to analyze the geometry of rough interfaces. The 
main result of this work is the dynamical roughness di- 
agram presented in Sec. IIVI based on scaling arguments. 
Section |V] gives the details of the performed numerical 
simulations that will be presented in Sec. IVII and give 
support to the proposed dynamical roughness diagram. 
Then, in Sec. I VIII we will present a discussion of the RM- 
RP crossover, the relation between the elastic string and 
particle chain models, the extensions of the roughness di- 
agram to the creep regime, and the implications of our 
results to numerical simulations with periodic boundary 
conditions. Finally, Sec. I Villi presents the conclusions of 
the present work. 



II. RANDOM PERIODIC SYSTEMS WITH 
LOCALIZED CORRELATION PEAKS 

As a model for a random periodic system with well 
separated length scales rf and M we focus our study on 
directed elastic interfaces described by a one component 
displacement field ii(r, t) with internal dimension d, r £ 
3?'', which satisfy an overdamped equation of motion 

jdtu{r,t)^cV^u{r,t)+Fp{u,r) + f + i^{r,t), (1) 

where 7 is the friction coefficient, c the elastic constant, 
and the uniform external force is given by F. The ther- 
mal fluctuations satisfy 

{V{r,t)) - 0, (2) 
(?7(r,t>?(r',t')) = 2jT5{t~t')S{r-r'), (3) 

and therefore the system asymptotically relax to the 
canonical thermal equilibrium at temperature T in the 
absence of the driving force /. The pinning forces are 
characterized by sample-to-sample fluctuations given by 

Fpiu, r)Fp{u', r') = A{u - u')S{r - r'). (4) 

In this paper we consider the case when the correla- 
tor function A(u) is a periodic function with correla- 
tion peaks localized in a range <C M at the values 
u = pM with M the periodicity and p any integer. We 
do not make a distinction between the so-called random- 
bond where we must enforce / du A(m) = and the 
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FIG. 1: (Color online) Schematic pinning force correlator 
for a random periodic system, with periodicity Ad and lo- 
calized correlation peaks with range rf. RB and RF stand 
for random-bond and random-field correlations (see the text 
for details). 

random-field cases (see e.g. Ref. [s^ since we are inter- 
ested in the depinning transition where these two dif- 
ferent static universality classes merge into a single one 
(in the creep regime, for f < fc and T > the distinc- 
tion must be done however, since the static properties 
can affect the intermediate length-scale physics^^*^) . In 
Fig. [T]wc schematically represent the shape of A(m). A 
periodic pinning force with fluctuations given by Eq. ^ 
arises naturally in numerical simulations of interfaces in 
random environments, when analyzing a system of trans- 
verse size M with periodic boundary conditions. We will 
exploit this fact to get most of our numerical results. 

Interestingly, as we will show later, we find that the 
roughness scaling for elastic lines in a random-periodic 
two dimensional potential with periodicity M in the di- 
rection of displacement also describe a chain of elasti- 
cally coupled particles in a one-dimensional non-periodic 
random potential, with the lattice spacing given by 
M. We argue that this connection is general, between 
a d-dimensional elastic manifold in a random-periodic 
medium and a periodic chain of {d — l)-dimcnsional cou- 
pled manifolds. The dynamics of a discrete chain of d- 
dimensional coupled manifolds can be described by 

^dtUn{r,t) = c[u„+i(r,t) -f- ■ti„_i(r, t) - 2u„(r,t)] 
+c\7^Un{r,t) + Gp{nM + u„,r) + F + i]nir,t), (5) 

where m„ describes the displacements of each manifold 
around the perfect position nM in the chain, c is a com- 
pression elastic constant, and G{u., r) an uncorrelated 
pinning force which is the same, independently of n, with 
a non-periodic short-range correlator of range rf. In this 
case, the correlations of the thermal noise are given by 
(77„(r,i)7?™(r,<)) = 27T<5„™5(r - r')S{t - t'). The con- 
nection between the physics of Eq. ^ for rf-dimensional 
manifolds and Eq. ([5]) for d — 1-dimensional manifolds 
is subtle since it involves a non-trivial coarse-graining in 
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TABLE I: Random-Manifold (KM) and Random-Periodic 
(RP) diaracteristic rougliness exponents for the three refer- 
ence states: equilibrium (EQ), depinning (dep), and fast-ffow 
(FF). In the static equilibrium case only the Random-Bond 
(RB) class is quoted. 



d= 1 


RM 


RP 


EQ 
dep 
FF 


Ceq = 2/3(RB) 
C = 1.25 
Cff = 1/2 


Ci^ - 1/2 

C''^ = 3/2 
= 3/2 



the direction of the periodicity which can produce extra 
terms in the equation of niotion i^^i^^i^" This issue will be 
discussed in more detail later. It is however plausible at 
this point that the resulting pinning force would display, 
i/ distortions are locally smooth, i.e. — Un\ <C M, 

well developed periodic correlations with a period M in 
the direction of the chain displacement with correlation 
peaks localized in a range rf <^ M , as the ones schemat- 
ically shown in Fig[TJ 



III. ROUGH GEOMETRY AROUND 
DEPINNING 



We focus our study on the geometrical observables that 
can be defined for the sliding manifold. The structure 
factor Sq is a very convenient quantity to study the ge- 
ometry of the manifold at different length scales and to 
locate the different crossovers] "^ We define it 
as 



e*«^M(r, t)dr 



(6) 



where we have chosen the particular direction x to be any 
component of r if the interface is governed by Eq. ([T]) 
and the direction of displacement (i.e. the direction of 
the chain) if it is described by Eq. ^ . 

The driven steady-state geometry at low temperatures 
is governed by three reference state a^^i^° : the / = equi- 
librium state, the f = fc and T = depinning critical 
state, and the fast-flow state / cxd. The particularity 
of these states is that above a microscopic length they 
have different self-affine geometries, i.e. the structure 
factor behaves as 



(7) 



where the power-law behavior reflects the lack of a char- 
acteristic length-scale in these states and C is the char- 
acteristic roughness exponent. The roughness exponents 
of the reference states are Ceq, C and Cff, respectively. 
These exponents can take different values in different uni- 
versality classes. While Ceq is different for the RB, RF 
and RP universality classes, C and Cff remain the same 
for RB and RF classes and they change for the RP class. 
Since we arc particularly interested in distinguishing the 



depinning and fast-flow roughness exponents of the RP 
class and the RM class, we use a superindex "RP" in all 
the exponent to indicate when the exponents belong to 
the RP case, and omit the superindex for the RM class 
(see Table HI). 

Furthermore, two important characteristic roughness 
exponents are the Larkin exponent Cl and the thermal 
exponent Cth , which are in general expected to appear at 
very small length-scales. The Larkin exponent is simply 
obtained by doing a first order perturbation expansion 
in the disorder, thus replacing it by a random uncorre- 
lated force. It yields Cl = (4 — d)/2 for lengths smaller 
than the Larkin length Ic, above which the naive per- 
turbation theory fails due to metastability. The thermal 
roughness exponent Cth = (2 — (i)/2 is defined as the one 
that appears in absence of disorder at finite temperature, 
and can be obtained exactly from the Edwards- Wilkinson 
equation. Interestingly, we will show later that both, 
Cl and Cth reappear at large length scales in the dy- 
namics of a RP system with localized correlation peaks: 

Cff — Cth, C — Cff — Cl. 

The steady-state geometry at small velocities can in 
general be described by velocity and temperature de- 
pendent crossover lengths separating different regimes 
of roughness. The corresponding roughness exponents 
arc however universal, velocity and temperature inde- 
pendent, and coincide with one of the aforementioned 
exponents. For velocities just above the RM depinning 
transition we have 




-(rf+2C) 

-(£i+2CFF) 



for q > 1/C 
for q < 1/C, 



(8) 



allowing to define the characteristic length C. For small 
velocities and vanishing temperatures C can be iden- 
tified with a velocity dependent divergent correlation 
length C ~ v~^/^^ . At / — )■ /+ and zero temperature 
we have C ~ (./ - ./c)"''r^ with w ~ (/ - fcf, and at 
/ = fc and small temperatures we have ^ ~ rp^u/p 
with ip a thermal rounding exponent such that v ^ T"^ . 
Since in this case Sq is governed by a single crossover 
length C we can write the scaling form 

d+2C 



(9) 



where the scaling function s{x) behaves as s{x) ^ 
^-{d+2CFF) for a; < 1 and s{x) for a; » 1. By 

plotting ^ YS V we can obtain a "geometrical roughness 
diagram" showing sectors with different roughness expo- 
nents at different observation length-scales I: C for I < C 
and Cff for I > C. Physically, C divides the small length- 
scales which arc dominated by the critical configuration, 
i.e. the unique v — steady-state solution of the equa- 
tion of motion for f — fc, from the large length-scales, 
which are governed by an effective Edwards- Wilkinson 
equation with a velocity dependent effective temperature 
dynamically induced by the disorder The physical ori- 
gin of this crossover is due to the fact that at small but 
finite velocity the renormalized disorder becomes a weak 
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perturbation at large enough length scales, acting effec- 
tively as a thermal-like noise in an Edwards- Wilkinson 
equation, with a short-range correlation time of order 
Tfjv and effective strength or "temperature" A(0)/w. 

In a random-periodic system with localized correlation 
peaks, as the ones described in the previous section, the 
scaling of Eq. © must be corrected to take into account 
the existence of the additional characteristic distance M. 
As we show in the next section, M induces new geomet- 
rical crossovers at depinning, separating the geometrical 
roughness diagram in more than two sectors. 



IV. DEPINNING ROUGHNESS DIAGRAM 
AND SCALING ARGUMENTS 

In this section we summarize our most important phys- 
ical results about the steady-state geometry of driven 
random-periodic systems with localized periodic correla- 
tion peaks. We present the depinning roughness diagram 
and heuristic scaling arguments describing the different 
crossovers. These arguments are corroborated numeri- 
cally and analytically in the following sections. 

In Fig. [2] we schematically show the geometric rough- 
ness diagram we find, by analyzing the structure factor, 
for a random periodic system with localized correlation 
peaks. It presents three roughness sectors, characterized 
by the roughness exponents of the RM depinning ^, the 
RM fast-flow Cff and the RP fast- flow Cff • Interest- 
ingly, unlike the RM case, for the RP system = (pp , 
and therefore there is no signature, in the steady-state 
structure factor, of the divergent length-scale expected 
for the depinning transition of a pure RP system. We 
discuss this issue later. Below a characteristic velocity 
vp the system crosses over, at a characteristic velocity- 
independent length Lp, from a small length-scale regime 
with a roughness exponent C, corresponding to the geom- 
etry of the RM critical configuration, towards a regime 
with an exponent C^^i corresponding to the RP criti- 
cal configuration. Above vp there are two crossovers, 
at the characteristic velocity-dependent length-scales ^ 
and Lm- The first crossover is from a regime charac- 
terized by the RM critical depinning exponent ^ to a 
regime with the fast- flow RM exponent Cff- The second 
crossover, observed by further increasing the observation 
length-scale, is from the RM fast-fiow regime to the RP 
fast-flow regime, the latter characterized by the exponent 
dominating the largest length-scales. 

The different length scales and roughness exponents 
shown in Fig. [2] can be obtained by analyzing the struc- 
ture factor. To illustrate how we obtain the roughness di- 
agram, in Fig. 131 a) we show a typical averaged structure 
factor for v > vp for an interface in a random periodic 
disorder medium with period M. Increasing the obser- 
vation length-scale (decreasing the wave vector q) we see 
different crossovers between different roughness regimes 
at the characteristic length-scales C and Lm- In Fig.[3l^b) 
we also show a typical structure factor for an clastic chain 
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FIG. 2: (Color online) Schematic roughness diagram at the 
depinning transition of a random-periodic system with period- 
icity M and localized periodic correlation peaks as a function 
of its steady-state velocity v. For v < vp, the geometry is of 
the RM (or domain- wall) class at small length scales I < Lp, 
while it is of the RP (or charge-density-wave) class at large 
length scales I > Lp, where Lp = Lp{M) but independent of 
th velocity v. For v > vp, the geometry is of the RM class 
at small length scales I < ^, where ^ = ^{v) is the RM depin- 
ning correlation length. By further increasing the observation 
length-scale a crossover between RM and RP fast-flow regimes 
of roughness occurs at the length Lm = Lm{M, v). Note that 
the large scale geometry is described by only one roughness 
exponent since fast-flow and depinning exponents coincide in 
the RP class, = Cff- The proposed scaling with the 

velocity i; and the periodicity M of the dynamical crossover 
lengths C('y)) Lp{M), Lm{v,M) (see text) is corroborated by 
analyzing the structure factor obtained from numerical simu- 
lations. 



with lattice spacing M moving in a one-dimensional dis- 
ordered medium, displaying identical regimes of rough- 
ness. This supports the argued connection between the 
geometrical properties of periodic chains of manifolds 
with internal dimensions d'm d-\- 1-dimensional pure ran- 
dom media and single interfaces of internal dimension 
d + 1 in d + 2-dimensional random-periodic media. 

The geometric roughness diagram for a RP system 
with localized correlation peaks is richer than for the 
RM system, which only displays one critical character- 
istic length-scale ^ above fc- This is due to the fact that 
M is an extra characteristic length in the problem, differ- 
ent from Tf, tmlike what occurs in CDW systems. Inter- 
estingly, the RP system we study thus contains the RM 
depinning diagram, with its C and Cff roughness sectors, 
for velocities larger than vp and lengths below Lm- In- 



deed, wc find C 



(independent of the periodicity 



M) which coincides with the divergent length scale of the 
RM depinning, being v and /3 RM critical exponents. We 
can thus say that below Lm for v > vp and below Lp 
for V < Vp periodicity effects, both spatial and temporal, 
arc not important. In other words, short lengths scales 
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FIG. 3: (Color online) Typical structure factor for (a) a mov- 
ing elastic string in a periodic random medium with peri- 
odicity M, and (b) an elastic chain with lattice spacing M 
moving on a non-periodic disorder medium. The velocity of 
these systems is larger than their characteristic velocity vp 
(see Fig. [21). The two systems display a critical RM depinning 
roughness at small length scales which crosses over to a RM 
fast-flow roughness at intermediate length scales, and then to 
a RP fast-flow roughness at the largest length scales. These 
regimes are separated by two characteristic dynamical length 
scales: the correlation length ^ and the periodicity-induced 
length scale Lm. 



sector for a CDW. 

To estimate Lm for v > vp we must be more care- 
ful. Indeed, the static argument given above of match- 
ing the width w of the shding manifold with M is in- 
correct in this case. This argument would give M ~ 
rf{i/Q^{LM/i)^^^ or Lm ~ i[{M /r f){lj , and 

by using ^ v'"/'^ we finally get Lm ~ ■y-'^ti-C/frF)/'^. 
Since in general > (^pF we get the incorrect result 
that Lm grows with the velocity, inconsistent with our 
data. The error in making such an argument comes from 
the fact that for v > vp the roughness of the interface 
is determined, above by the temporal correlations 
of the pinning force (when seen from the moving inter- 
face). Indeed, despite the fact that w < M and that 
the rcnormalizcd disorder is already weak at peri- 
odicity effects are relevant above a certain length beyond 
which the manifold has not time to relax all its modes 
after moving by a distance M. The steady-state geom- 
etry of the moving system thus probes the periodic cor- 
relations of the pinning force. We must thus compare 
the typical relaxation time in the RM fast-flow regime 
^ for C < ^ < Lm, with the "time 

of flight" tm = M/v, being zff = zth = 2 the dynami- 
cal exponent of the fast-flow RM class, z the dynamical 
exponent of the critical RM regime, and Tc a microscopic 
time. If these times equate at Lm, 
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ZFF 



(11) 



(12) 



This result with x > describes well our numerical data, 
as wc show later. 

Having L]\i and Lp wc can now determine the charac- 
teristic velocity vp of the roughness diagram, defined as 
Lm{vp) = Lp. We get, 



decorrclate rapidly, spatially and temporally, and can not 
sense the periodic correlations of the pinning force. 

To understand Lp we can make the simple scaling hy- 
pothesis that periodicity effects start to be important 
when the average transverse size or width of the critical 
manifold is of the order of the periodicity M, w ^ M. 
Since at small lengths I, the width grows as w = rf{l/lc)'' 
(both for V < Vp and v > vp) we can compute the char- 
acteristic length Lp by stating that M ~ rf{Lp/lc)''- We 
thus obtain that, 



Lt 



(10) 



independent of the velocity. This describes well our data 
as a function of M,v. Physically, Lp thus represents 
the length at which the critical configuration starts to 
see the periodic spatial correlations of the pinning force. 
Note that for the CDW case we have M ^ rj and thus 
Lp ~ Zc- Therefore we would not observe this RM critical 



and therefore. 



Vp 



(13) 



(14) 



It is worth noting here that while vp decreases, Lm 
and Lp increase with increasing M . This means that for 
large enough M the RM sector of the roughness diagram 
of Fig. [5] grows and in practice the system behaves as 
a RM system. Conversely, for small M the RP sector 
grows and dominates the behavior at small velocities. 

We also note that below up and above Lp we expect 
to observe RP or CDW-like depinning, with a divergent 
correlation length ^P ~ {f — fc)^'^ ■ However, unlike the 
RM case, the divergent length does not manifest itself as 
a crossover between roughness regimes of the structure 
factor, since = Cff • This is consistent with the fact 
that the roughness exponent = Cl = (4 — (i)/2 for the 
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RP appears in FRG calculations from the generation of 
a random-force in the renormalized pinning correlatorj^S 
In other words, the pinning forces acting on pieces of 
size Lp are essentially uncorrelated, and the model thus 
effectively becomes the Larkin model, with a roughness 
exponent = (4 — d)/2. In this respect, in Sec. IVII we 
show that Cpp = from numerical simulations. 

The roughness diagram of Fig. [2] appears to be valid at 
small finite temperatures within the "thermal rounding" 
regime, as we find numerically. In this regime the effect of 
the temperature translates into a finite velocity v ^ T"^ 
at f = fc but does not affect the large-scale roughness 
regimes. The depinning roughness diagram of Fig.[5]thus 
remains the same, whether the velocity is originated by 
driving force, small temperature or both. 

In the following sections we describe our numerical sim- 
ulation method and results supporting the roughness di- 
agram of Fig. [3] and the scaling for the different crossover 
lines. 

V. DETAILS OF NUMERICAL SIMULATIONS 

We present here a detailed description of the numerical 
methods we use to study the RP system with localized 
correlation peaks. For simplicity we analyze low dimen- 
sional manifolds but our results remain qualitatively the 
same for higher dimensions. The d = 1 case turns out 
to be on the other hand the most stringent case for our 
general arguments. 

We study the motion of an elastic string in a disordered 
environment described by Eq. ([T]) in d = 1 (D = 2). In 
order to numerically solve Eq. ([T]) for the elastic string we 
discretize the D = 2 embedding medium in the longitu- 
dinal z-direction in L segments of unit size, keeping the 
transverse displacement field u(z) as a continuous vari- 
able in the x-direction. The discrete system of equations 
read 

^dtu{z,t) = c + 1, <) + it(z — 1, t) — 2it(z, i)] 

+ Fp{u,z) + f + ri{z,t), (15) 

with z an integer. Periodic boundary conditions of size 
M (resp. L) are imposed for the transverse (resp. lon- 
gitudinal) system sizes. Besides avoiding boundary ef- 
fects, this model presents several advantages which have 
been exploited in various ways for non-periodic sys- 
tems The critical force and critical config- 
uration for such finite systems can be determined for 
each sample in polynomial time with arbitrary precision 
by exploiting the Middleton theorcmsj^ Moreover, the 
complete sequence of metastable states below threshold, 
and in particular the one dominating the creep motion 
at low temperatures can be determined for each particu- 
lar sample exactly, by generalized, Middleton-like, theo- 
rems 4^1^ This method thus allows for a well controlled 
analysis of key properties such as the critical force statis- 
tics, the critical exponents of the depinning transition, 
and the different roughness crossovers. 



To show the generality of our results we also study for 
comparison the problem of an elastic chain by solving 
Eq. © in D = 1 and d = 0, 



7i9tu„(t) = c[u„+i(t) 4- u„-i(t) - 2M„(t)] 

+ Gp{nM + Ur,)+ f + Tln{t). (16) 



As mentioned above, Fig.[3]shows that the elastic string, 
described by Eq. and the elastic chain, described by 
Eq. (|16p display the same roughness crossovers around 
depinning. We argue that this geometrical equivalence 
is general, between the d dimensional manifold in a 
D = d + 1 disordered medium with period M and the 
periodic chain oi d — 1 dimensional elastically coupled 
manifolds with lattice spacing M in a Z) = d dimen- 
sional disorder medium. Therefore we study in details 
the case of the elastic string and translate appropriately 
our results to both kinds of systems in any dimension d. 
The resulting discrete system of equations for the elastic 
chain in D = 1, Eq. (jl6p are indeed similar to the ones 
for the string in 13 = 2, by identifying the discrete values 
of z for the particles of the string with the index n for 
the particles of the chain. The main difference between 
the two systems are the pinning force correlations. While 
the pinning force on the string is uncorrelated for differ- 
ent values of the labelling variable z, the pinning force on 
the chain is correlated for different values of the labelling 
variable n, since in the latter case the particles visit the 
same disorder as they move. This difference can be bet- 
ter appreciated by remarking that the equations for the 
elastic chain are equivalent (by interpreting n as z) to the 
ones of a tilted elastic string in a Z? = 2 medium with 
columnar disorder, being 9 = ta.n^^{M/L) the imposed 
tilting angle. The result of Figs. |3] is thus non-trivial and 
suggests that the roughness diagram of Fig. [2] is general 
for elastic system with pinning forces displaying localized 
disorder correlation peaks. 

The equations of motion, Eqs. and (|16p . are in- 
tegrated using Euler method with a time step 5t = 
0.01. We set 7 = 1, c = c = 1, r/ = 1, and 
a disorder strength A(0) = 1. A different choice of 
this microscopic parameters does not qualitatively alter 
our results. The continuous random potential for the 
string V(u,z) ^ — J du Fp{u,z) is modelled by L cubic 
splines passing through M regularly spaced uncorrelated 
Gaussian numbers points. For the chain, the potential 
V{u) = — J du Gp{u) is numerically generated with ran- 
dom spline passing through L x M regularly spaced un- 
correlated Gaussian numbers points. Disorder average is 
done by averaging over different realization of the gaus- 
sian random points. Using these disorder potential mod- 
els, when M ^ Tf the corresponding pinning forces dis- 
play periodic correlations with localized peaks in a range 
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FIG. 4: (Color online) Scaling of the structure factor for dif- FIG. 5: (Color online) Scaling of the structure factor for dif- 
ferent M values at T = and / = 1.2/c. The longitudinal ferent M values ai f — fc and T = 0.5. The longitudinal 
system size is L = 512. (a) Raw data, (b) Scaled data. system size is L = 512. (a) Raw data, (b) Scaled data. 



VI. NUMERICAL RESULTS 

In this section we show and discuss the numerical re- 
sults for the characteristics lengths, roughness exponents 
and characteristic velocities, appearing in the geometri- 
cal roughness diagram of Fig. [21 We describe separately 
the different crossovers for v > vp and v < vp. 



A. Roughness Crossovers for v > vp 

We start by discussing the two crossovers observed in 
Figs. IHa) and ^h) for v > vp, at the characteristic 
lengths ^ and L m respectively. One observes three rough- 
ness regimes for given values of L and M. They corre- 
spond to the three regimes observed in Fig. [5] for veloci- 
ties above vp. Increasing the length-scale (decreasing the 
wave vector q) the local roughness exponent changes from 
C w 1.25, to Cff 0.5, and finally to Cpi? = Cl = 1-5. 
The first two roughness exponents are characteristic of 
the RM depinning, and we can identify the crossover 
length ^ with the divergent correlation length of the RM 
depinning down to vp. The second crossover length Lm 
is proper to our system, separating the fast-fiow regime 
of the RM class from the one of the RP class. 

For pure RM models the value ( w 1.25 was obtained 



numerically before by very different methods ;^^'^ and in 
particular by exact algorithms i^i*^ Two loop renormal- 
ization group calculations for the RM class are required 
to get values that are consistent with this result 3^ The 
exponent ^ff was obtained by numerical simulations)^ 
and by analytical argumentsi^ As described above the 
physical meaning of the appearance of (^ff is that at 
large length-scales the velocity v becomes very important 
and disorder effectively acts as an spatially uncorrelated 
time-dependent perturbation with short-range temporal 
correlations in a range rj/v. The strength or "effective 
temperature" of this effective noise is thus proportional 
to A(0)/'i;i^ Equation ^ then effectively becomes an 
Edwards- Wilkinson equation for which it is straightfor- 
ward to show that the steady-state roughness exponent 
is Cth = (2 - d)/2 for d < 2, and Cth = for d > 2. 
For d = 1 we have Cff = Cth = 1/2, consistent with our 
numerical result for the present system. 

As described in Sec. IIVI the crossover at Lm for v > 
f p occurs when the relaxation time of the string in the 
RM fast-flow regime becomes of the order of the time of 
flight M/v. The resulting scaling of Lm with M and v 
involves several exponents of the RM class, and hence it is 
a good test for the validity of our scaling arguments. For 
fixed external force and temperature and changing the 
transverse size M only the crossover around Lm changes 
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FIG. 6: (Color online) Scaling of the structure factor for dif- 
ferent f > fc values at T = and M < L'' . The longitudinal 
system size is L = M = 512. (a) Raw data, (b) scaled data 
around ^, and (c) scaled data around Lm. 



FIG. 7: (Color online) Scaling of the structure factor for dif- 
ferent T > values at f = fc and M < Lf' . The longitudinal 
system size is L = M — 512. (a) Raw data, (b) scaled data 
around ^, and (c) scaled data around Lm. 



in the structure factor, as observed in Figs.SJJa) andlSJa). 
Thus, for the crossover at large length scales the structure 
factor can be written, for g <C as, 

^,L;f = G{qLM) (17) 

with G{x) ~ x'^^+^^ff) for a; < 1 and G{x) ^ 
^-(i+2Cff) foj. j; > 1. Using Eq. dTT]) for Lm, and since 
the velocity is fixed by / and T, we get the following 
scaling formula, 

S',M-(i+2«'''')/^ff ~ G{qM^/'""']) (18) 

In Figs, mb) and [Hb) we test this scaling prediction 



as a function of M for a fixed velocity by plotting 
S',M-(i+2frF)/^FF vs gAfi/'^ff for different values of M. 
In Fig.mthe velocity is produced by a force above thresh- 
old at zero temperature while in Fig. [5] the velocity is 
produced by a finite but small temperature at / = /c. 
In both cases we find that the scaling form proposed col- 
lapses well the different curves by using z-pp = 2 which 
corresponds to the dynamical exponent of a RM model 
at large velocities (i.e. the dynamical exponent of the 
Edwards- Wilkinson equation with the disorder-induced 
Langevin-like noise). As expected, deviations from the 
good collapse are observed only at large q, where the 
presence of the extra characteristic length ^ invalidates 
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the simple scaling of Eq. . We also note that in the 
non-scaled data in Figs. |H[a) and \5la) Sq becomes M- 
independent for q ^ This is consistent with the 

fact that ^ does not depend on M but only on the veloc- 
ity, ^ ^ v~^/^ near depinning, unlike Lm which depends 
on both, V and A/. 

In order to study the velocity dependence of the struc- 
ture factor and its crossover Icnghts we have applied both 
different driving forces / > /c at T = and small tem- 
peratures for / = fc- In Figs.|ni[a) and [71[a) we show Sq 
as a function of the force and temperature, respectively. 
The crossover around ^ can be described, for q ^ ' 
with the scaling relation 



(19) 



with the function G{x) ^ x (i+2Cff) for a; ^ 1 and 
G{x) ^ x"'^^^'''' for x ^ 1. This scaling form depends 
on force and temperature only through v. To get explic- 
itly these dependencies we can use, for small v and f > fc 
that v'^ p for f> fc, with f ={f- fc)/fc the reduced 
force, and v - T"^ for f = fc and small TM Figures E^b) 
and [7{b) show the respective scaling forms around ^. 
However, this scaling form is valid up to the scale Lm 
where periodicity effects are important and G{x) is not 
longer universal. The crossover of G{x) to RP fast-flow 
can be written as 



G{x)r.x-^'^^^^^^G(—\ 



\XM J 



(20) 
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FIG. 8: (Color online) Scaling of the structure factor for dif- 
ferent M values in the v < vp regime. Data correspond to 
T = 0.001 and L = 512. (a) Raw data, (b) Scaled data. 



where Xm — v ILm and the new function G{y) ^ 
y-{\+K^l) for y <C 1 and G(y) ~ y-(i+2CFF) foj. y > 1. 
Using Eq. ([TT|) one can write, for fixed M and for g <C 1 
that the structure factor behaves as 



with 



v-'^Giqv-^), (21) 



K = 2(C-CFF)i^//3 + (l + 2CFF)x- (22) 



This form describes the crossover between RM and RP 
fast-flow regimes of the structure factor. Figures[6Kc) and 
[71Jc) show this scaling form when the velocity is generated 
by finite drive at zero temperature or by a small finite 
temperature at f — fc- 

Therefore, in Figs. IH] and [7] we test the latter scaling 
prediction as a function of / for f > fc and T = 0, and 
as a function of T for / = /c by using the known values 
/3 = l/3^i^ V = 0.154^ Since Sq has two characteris- 
tic lengths Lm and ^, we show separately the collapse 
around the two crossovers. In Figs. IHJb) and [Zljb) we 
show the collapse around ^ and in Figs.lHJc) and ^c) the 
collapse around Lm, for the same set of curves Sq (/, T) of 
Figs.jBJa) and [DJa), respectively. The collapse obtained 
by using the known values of ip, (3 C,ff, z and zff fully 
supports our interpretation of the two crossovers. 



B. Roughness crossovers for v < vp 

As shown in Fig. [2l at vp we have L]\j{vp) = £,{vp) = 
Lp, where Lp is independent of the velocity. The RM 
fast- flow regime thus disappears at vp . This implies that 
periodicity effects are already generated by static pinned 
configurations, and not dynamically as described in the 
derivation of Li\j. Lp is precisely the length at which, 
for V < Vp, the width of a typical RM critically pinned 
configuration reaches M. Below Lp, periodicity effects 
are absent, and the typical critical configuration is not 
sensitive to M. Just above Lp, the critical configuration 
crosses over to the RP class. To prove the existence of 
this crossover we show in Fig. |8] that the structure factor 
of a string with L > Lp and v < vp can be written as 



SqM-^'+^<^/i = H{qM'/i 



(23) 



with H{x) - for X < 1 and H{x) - ^-(i+^C) 

for a; 1. We have used L ~ Af^^^ and the previously 
known values of and (. The collapse of Fig. [5] thus 
supports our interpretation of this crossover. 

As discussed in Sec.lIVIand shown in Fig. [51 we expect 
a second crossover at a length above Lp, represent- 
ing the dynamical correlation of the RP depinning. The 
regime between Lp and thus represents the RP critical 
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regime, with the roughness exponent of the RP critical 



configurations. As discussed in Sec. llVi C 



c, 



RP 
FF ' 



and 



due to this the RP depinning correlation length can 
not be detected by analyzing Sq. This suggests that sub- 
tle geometrical measures are probably needed to locate 
Cp^ 



VII. DISCUSSION 



proportional to A{0)/'yv. Within this approximation, the 
equation of motion becomes linear, and its solution in 
Fourier space is, for a particular component q = q.z of 
the wave vector q, 

Jo 

(25) 

The instantaneous structure factor is thus given by 



In this section we discuss, by analytical arguments and 
additional numerical simulations, the results obtained in 
the previous sections. We first discuss, in section IVII Al 
how to analytically calculate the crossover from random- 
manifold to random-periodic fast flow regimes of rough- 
ness at the length-scale Lm, which is absent in the pure 
Random-Manifold or Random-Periodic depinning. The 
numerical results of Sec lVII support the depinning rough- 
ness diagram of Fig. [5] which is found to be the same 
for interface pinning potentials with periodic correla- 
tions and for periodic elastic systems such as chains. In 
Sec. IVII Bj we discuss this interesting equivalence between 
the geometry of periodic elastic systems in random pin- 
ning potentials and interfaces in random-periodic pinning 
potentials. We give analytical and numerical arguments. 
In Sec. IVII Cl we discuss how to extend the roughness dia- 
gram of fig. [5] to the creep regime in the low temperature 
limit, based on the effects of periodicity in the statics 
and depinning. Finally, in Sec. IVII^ we discuss the im- 
plications of the roughness phase diagram for numerical 
simulations of elastic strings with periodic boundary con- 
ditions and the thermodynamic limit. 



A. Random-manifold to Random-periodic 
crossover in the fast flow regime 

We show here how the crossover from random-manifold 
to random-periodic fast-flow and their respective rough- 
ness exponents can be analytically computed. We con- 
sider here a simplified model of disorder perturbativcly in 
the large velocity and small temperature limit of Eq. ([T51) . 
We will show this approach yields correct results for 
the velocity and periodicity dependence of the crossover 
length, which can be therefore identified with Lm in 
Fig. [5J when using the renormalized friction, disorder 
strength and temperature at the length-scale ^. 

At high velocities and small temperatures, we approx- 
imate the pinning force Fp{u,r) by Fp{vt,Y). In doing 
so, we assume, a priori, that the disorder-induced and 
thermally-induced displacements are small in the regime 
we arc interested in. In this approximation the pinning 
force becomes an effective tlicrmal-likc noise fj{t, r) = 
Fp{vt,r) with temporal correlations given by 



fi{t, r)?7(t', r') = A{v{t - t'))5{r - r'). 



(24) 



Since A(x) has a spatial range rj, the range of temporal 
correlations of rj is rf/v and its effective temperature 



S,{t) = {\u,{m = Sf'^it) + 5f (t), (26) 

where the first contribution is the Edwards- Wilkinson or 
purely thermal structure factor 



S^^it) = ^ (l - e--V.) , 



(27) 



and the second contribution comes from the disorder- 
induced noise ry. 



(28) 

To proceed we assume a particular periodic correlator 
A(u) = A{u + nM), with n an integer. A simple choice 
for A{u) having sharply localized peaks at m = nM is 



JO 



A(u) = Ao ^ S{u - nM), 



(29) 



This kind of disorder arises physically from a random 
distribution of identical point like pinning centers acting 
on a very thin interface, in the limit rf with the 

constraint jj^^ dx A{x) = Aq and Aq a constant mea- 
suring the strength of the disorder. With such disorder 
the fast-fiow contribution to the structure factor can be 
easily integrated to get 



An 



V 7 1- e-^' 



(30) 



Then, the total instantaneous structure factor can be ex- 
pressed as 



^,(0 = 5^(0(1+ 



1 



2^VT 1 - ^-Icq^Mhv 



(31) 



Since we are interested in the steady-state we take the 
t ^ oo limit to obtain the steady-state structure factor, 



s ^ (ti ^ 

cq^y 2-fv 1 - e-2(9;A/)^ 



(32) 



which presents the characteristic length Im = ^JMcfjv 
or characteristic time tm = M/v. 

For large length-scales such that q <C l^} we have 



T 



■}q ~ —q 

c 



Aq _4 



(33) 
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where the velocity v does not appear explicitly. At T = 
0, Sq ^ q~^, implying that the large scale roughness is 
identical to one of the Larkin model Cl = (4 — (i)/2 if 
d < A and Cl = otherwise. This is consistent with 
our finding Cff = Cl = 3/2 for the d = 1 case (string) 
in a random periodic medium and also for the elastic 
chain. If T > we could have a roughness crossover at 

qr = \J 4^^, where the two terms in Eq. ([55)) become 
equal. However, for this roughness crossover between Cth 
(with Cth = (2 — (i)/2 for < 2 and Cth = otherwise) 
to Cl to be observable we must require qr ^ or 
equivalently T ^ Aqc/v. Since Aq/v can be identified 
with the shaking temperature Tgh of Ref . [69l the condition 
for observing this crossover in this regime reads Tshiv) ^ 
T, meaning that disorder-induced fluctuations must be 
much smaller than thermal fluctuations. 

U q ^ l^j on the other hand, we have 
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FIG. 9: (Color online) Ground state geometry of the elastic 
string (open symbols) and the one-dimensional chain (closed 
symbols) problems as obtained with transfer matrix calcu- 
lations. The data are represented in a scaled form and the 
closed symbols are shifted upwards for clarity. 



which is equivalent to the Edwards- Wilkinson structure 
factor at an effective temperature Tcs(v) = Tsh(v) + T, 
yielding a roughness exponent Cth- 

From ([Ml) and ([55[) we thus see that if the temperature 
is small compared to Tghiv) = Aq/2v there is a rough- 
ness crossover at the length-scale Im from Cff = Cth to 
Cff ~ Cl when increasing the observation length-scale. 
This is in agreement with the roughness diagram of Fig.j^] 
for V > vp and lengths above ^. By comparing Im and 
Lm we see that both quantities have the same M de- 
pendence, since zff = 2. The explicit velocity depen- 
dence, even if it is a power law in both cases, is differ- 
ent. In this model we find Im ^ (7''^)"^^^! instead of 
Lai i)^i/2+(''//3)(2/2-i) which describes well our data 
and was predicted in Eq. (jlip by pure scaling arguments. 
We also note in this respect that the structure factor at 
small q predicted by the model appears to be velocity 
independent and temperature dependent, in contrast to 
what we predict by scaling arguments and what the data 
presented in Sec. IVII supports. These differences can be 
directly attributed to the incorrect use of the bare fric- 
tion constant 7, disorder strength Aq and temperature T 
in our perturbation theory, instead of using their renor- 
malized velocity-dependent values 7(w), Ao{v) and T at 
the length scale ^. To prove this we first note that at 
the length-scale ^ we have 7(1; )w = (/ — fc) = 
and therefore 7(f) = Then, by replacing 7 

by 7(f) we get the same velocity and periodicity depen- 
dence for Im and Lm- This justifies the identification of 
the crossover predicted with the present model with the 
one from the RM to the RP fast-flow regimes observed 
in the simulations, and estimated in Sec. llVl bv physical 
arguments. 



B. Elastic String vs Elastic Chains 

We discuss here the equivalence observed between the 
geometry of an elastic line in a two-dimensional RP 
potential and the one of a periodic chain in a one- 
dimensional non-periodic random potential. We argue 
that this connection is general, between D-dimensional 
thin interfaces transversely displacing in random-periodic 
D-l-1 dimensional spaces and a periodic array of (D-1) di- 
mensions coupled interfaces in a D-dimensional random 
medium. The connection between the two systems is 
however not trivial, since there is no exact mapping be- 
tween these two systems. We describe first the case of the 
statics, which can be discussed in terms of the replicated 
Hamiltonian and complemented with additional transfer- 
matrix numerical calculations, and then the dynamics. 



1. Statics 

We start by analyzing first the static problem, i.e. 
F = and T = 0, of the elastic line in a two dimensional 
random-periodic disorder and then compare it with that 
of a one-dimensional elastic chain over a random poten- 
tial. 

Let us consider an elastic line described by the unival- 
ued function u(z), where w is a function in the transverse 
direction, and z gives the longitudinal direction. The 
elastic contribution to the Hamiltonian is 

H.^^Jdz id,u)\ (35) 

while the disorder contribution in term of the line density 
p{x, z) is 

Hdis — / dx dzV{x, z)p{x, z). (36) 
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Using the integral representation of the 5 function one 
can write the density as 

pix, z) = S{x - u{z)) = -L /" dA e*^(^-"(^» . (37) 

271 J 

In order to consider the periodicity of the system in 
the x-direction, the disorder potential can be written as 
a sum over periodic images as 



Vix,z) = J2vix-pM,z), 



(38) 



Each 'mass' has a finite internal width ui and the center of 
consecutive masses are separated a distance a. The posi- 



Ra- 



tions of the particles are given by Xj 
where x^ = ja is the nominal position in the unperturbed 
lattice and uj is the displacement. In order to treat the 
model one should go from the Uj variables to a continuum 
formulation. Through this relabeling process, the decom- 
position of the density in this one dimensional problem 
leads to a density field 



X — a;,- 



where V{x, z) is the disorder potential defined in the in- 
terval X G [—M/2, M/2]. Using this periodic potential 
one can write for the disorder Hamiltonian that (see the 
Appendix) 



Hdis — — 



A 



2TM 



(39) 



a.h K 



where K = 2Tin/M . This last expression is strictly iden- 
tical to the one of a periodic system for which K are 
the vectors of the reciprocal lattice, as we show in the 
following. 

In order to describe the one-dimensional periodic chain 
in a one-dimensional disordered potential, let us consider 
an elastic chain of average lattice space a, which corre- 
sponds to average density po = !/«• One can imagine a 
chain composed of masses and springs of constant length. 



Pa 



1 — dxu{x) + 



^iK[x — u{x)) 



, (40) 



where the last expression is valid for 9^:1* <C 1, and K = 
2Tm/a are the reciprocal lattice vectors. The continuum 
field 



u{x) 



2-K/a 



27r 



(41) 



is valid for a; ^ a. For details on the relabelling process 
see Ref. [gO, especially Appendix A. 

Considering the decomposition of the density and using 
a replica formalism, the replicated disorder Hamiltonian 
can be finally written as (see the Appendix) 



2T ^ 

ah 



dx 



dxu°'{x)dxv}'{x) 



■iK(u''(x)-u*(x)) 



(42) 



which, with the exception of the dxu"" {x)dxU^ [x) term, 
is formally equivalent to Eq. ([M]). However, the original 
one dimensional chain in a one dimensional disorder po- 
tential, as studied by Cule and Hwa,— contains also a 
term proportional to dxu{x). This term is irrelevant at 
large length scales in dimensions d > 2, and gives a finite 
shift of the correlations function in d < 2. Indeed it is 
commonly accepted that this term does not change the 
roughness exponent j ^^i^^ 

Then, when comparing the periodic disorder case to 
the periodic chain, the period of the disorder potential M 
becomes the average distance between neighboring parti- 
cles. Thus, when the fiuctuation of the particles becomes 
of the order of the average distance, the system enters 
a RP or CDW regime. Then, in the static limit, struc- 
tural fluctuations at large length scales are given by the 
roughness exponent associated to the static CDW prob- 
lem, = 1/2. At finite M > 1 a crossover appears at 



a given length scale, corresponding to the scale for which 
the disorder induced fluctuations become of the order of 
the periodic box. 

The geometrical equivalence argued above can be 
tested directly by transfer matrix calculations for the 
ground-state of a one-dimensional chain and the one of an 
elastic string. In Fig. [S] we compare the structure factor 
of both systems, for different periodicities M :§> r/. We 
see almost no difference, meaning that the extra terms in 
the replicated Hamiltonian for the chain are irrelevant. 
We also find that the structure factor has a crossover at 
a characteristic length ip = Lp{M) between a regime 
with the RM equilibrium roughness exponent i^eq = 2/3 
at small length scales, to a regime with the RP equi- 
librium roughness exponent Cgq = 1/2 at large length 
scales. In Fig. [9] we show that the structure factor is well 
described by the scaling formula 

5,^q-(i+2«-«)7?(gL0), (43) 
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with a crossover length Lp ~ M^/'-'^'S. The function 
H{x) is such that for small x it behaves as ^^'•^'^'^"^eq) 
The crossover at Lp can be understood in the same terms 
as for the crossover length Lp at depinning. Lp is in this 
case the length at which the width w{L) ^ L^'^'s of the 
interface at the ground-state becomes of the order of the 
periodicity i(j(Lp) ~ M. As for depinning when M 3> rf 
the structure of the system at equilibrium is identical to 
the RM one, and crosses over to the RP one at large 
length-scales. 



at intermediate and small length-scales for which the RP 
systems behave as RM (or RF) systems. Since at large 
length-scales the effect of periodicity should always ap- 
pear, it is plausible to connect the crossover lengths Lp 
at / ^ /c and Lp at / = by a dynamic crossover length 
lp = lp{f,M) in the creep regime. The width of the 
interface at Lp is thus given by 



if / < /o and 



M - (L^) 



0\Ceq 



(47) 



2. Dynamics 

One can do a similar comparison for the dynamics of 
the clastic line in the periodic disorder potential and the 
periodic chain in disordered medium. As usual, we model 
the motion of an elastic string in a disordered environ- 
ment by means of the overdamped equation of motion 
Eq. pS]) . The pinning force is 



(44) 

If now one use the Martin-Siggia-Rose formalism as in 
Ref. [tO, then it can be shown that Fp interacts only lin- 
early with the operator iizt- Thus one can replicate over 
different times in order to obtain an averaged expression 
similarly to what was done in the previous section, which 
results in a term proportional to 



where we used that 



Fp{x, z)Fp{x', z') = AFS{x - x')S{z - z'). 



(45) 



(46) 



Again, this is equivalent, up to a factor proportional to 
dxu{x), to the case of the one dimensional periodic sys- 
tem^ The main difference is the convective term, but it 
can be shown that it is irrelevant, since an arbitrary shift 
u{x) u{x)+f{x) leaves the disorder term unchangedi^ 



C. Roughness Diagram in the Creep Regime 

The roughness diagram of Fig. [5] for depinning and 
the results of Fig. |9] can be combined to infer a rough- 
ness diagram as a function of the driving force, includ- 
ing the expected crossovers in the creep regime, / < fc- 
As it was shown in Ref. [13 for pure RM (or RF) sys- 
tems, below the depinning threshold fc a characteristic 
length Lopt = Lopt(/) exists. This length grows with 
decreasing / separating the small length-scales described 
by the equilibrium (/ = 0) roughness from the depinning 
(/ = fc) roughness. Since M ^ rf, we can expect to ob- 
serve the same behavior for the random-periodic system 



M L 



Ceq 
opt 



Jopt 



(48) 



if /p < / < /c, with fp a new characteristic force, defined 
by the condition lp{f^,M) = L%{M). Therefore 



M 



i/C 



Jopt 



l: 



Ceq 
opt , 



(49) 



with /p < f < fc- In the diagram of Fig. [TO] we schemat- 
ically show the crossover length lp separating the RM 
from the RP depinning roughness in the creep regime. It 
shows several sectors, including the equilibrium, depin- 
ning and fast-flow geometries of both the RM and the 
RP (note however that the fast-flow and depinning RP 
regimes have the same roughness exponent). 

Note that at variance with Lm, lp can be obtained by a 
static argument similar to the one used for Lp since the 
velocity vanishes rapidly in the low temperature creep 
regime. The system has thus enough time to relax all 
its (non-zero) modes in the time of flight tm = M/v. 
Note in this respect that the largest length-scales obey an 
Edwards- Wilkinson equation with correlated noise, and 
thus their relaxation times are governed by a dynamic 
exponent, contrarily to the zero-mode displaccm 



D. Periodic boundary conditions and the 
thermodynamic limit 

We present here a discussion concerning numerical sim- 
ulations and the proper thermodynamic limit. It is well 
known that when using numerical simulations to describe 
interface depinning and creep, a main shortcoming of the 
numerical method is coming from the difficulty in taking 
the thermodynamic limit M ^ oo and L — >■ oo. In or- 
der to perform a consistent finite-size analysis one has to 
carefully specify how both M and L should tend to infin- 
ity, as different prescriptions for the aspect ratio scaling 
lead to very different results. 

It has been shown that the sample-to-sample fluctu- 
ations of the critical force drastically change with the 
ratio M / For very small M compared to L^ period- 
icity effects are important and the distribution of critical 
forces is Gaussian, while at very high values of M the 
critical force is dominated by extreme values and its dis- 
tribution becomes of the Gumbel form. In the first case 
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FIG. 10: (Color online) Depinning roughness diagram of 
Fig. [5] extended to the creep regime < / < /c, at small 
temperatures. The diagram displays roughness exponents of 
the RM and RP classes at equilibrium, depinning and fast- 
flow. 



values of the different roughness exponents. The geomet- 
rical roughness diagram of Fig. [5] shows clearly that this 
is indeed the case, and gives at the same time an answer 
to this problem. It shows that for a fixed value of M, the 
optimal aspect ratio scalings are L = Lp{Al) for v < vp, 
equivalent to those proposed in Ref. l63l . but a different 
aspect-ratio scaling, L = Lm{v, M), for v > vp. Figure [2] 
thus shows that using a velocity independent prescrip- 
tion L = Lp{M) ^ Af^/^, which works for the critical 
configuration, would always give inconsistent results at 
all non-zero velocities in the thermodynamic limit, since 
wp — )- when M ^ oo and then L ^ Lm at a fixed 
V. Therefore, by increasing M within this prescription 
the system would eventually display periodicity induced 
eifects at any finite w, inducing an artificial crossover as 
a function of the velocity in non-periodic systems. This 
crossover induced by periodicity is on the other hand 
physically interesting for RP systems with localized cor- 
relation peaks. 



VIII. CONCLUSIONS 



the mean critical force is always bounded while in the 
latter grows logarithmically with Ad and is thus infinite 
in the thermodynamic limit. For the aspect ratio scaling 
M ~ it was shown that the mean critical force is fi- 
nite and well defined in the thermodynamic limit where 
the system displays pure RM behavior. In this last case, 
the critical force distribution is between Gaussian and 
Gumbel. The aspect ratio scaling M = kL'' leaves how- 
ever open the question of what is the optimal value of 
k for avoiding effects induced by the transverse bound- 
ary conditions. This has motivated the use of a differ- 
ent method for calculating steady-state properties of the 
same system, in which the control parameter is not the 
force but the mean velocity of the manifold by replac- 
ing the driving force by a uniform spring with constant 
m^, / — ?► ■m'^{vt — u(r)), pulling the whole manifold at a 
constant speed v in a transversely infinite medium. This 
model allows for more direct comparisons with analyti- 
cal calculations;^"— and has the advantage, compared 
with the force-controlled model, that the characteristic 
length Lm induced by the parabolic moving potential is 
controlled only by the spring constant ~ 1/m, and 
not by the velocity-dependent geometry of the manifold. 
When modelling a system for which the spring has not a 
physical origin the correct scaling for the spring constant 
is simply L„i ^ L in this case and thus very small values 
of m are required in the large-scale limit. 

If one wants to stick to the constant-force model the 
most natural empirical choice for the aspect ratio scaling 
is 'w{L) ?a A/, where w is the average width of the mani- 
fold)^ The problem with such prescription is that w{L), 
being an integrated quantity w ~ j dq Sq, has a compli- 
cate dependence with the mean velocity of the manifold 
through the velocity-dependent correlations lengths sep- 
arating different regimes of roughness, and through the 



We have studied numerically the depinning transition 
of driven elastic interfaces in a random-periodic medium 
with localized periodic-correlation peaks in the direction 
of motion. We have obtained a dynamical roughness di- 
agram which contains, at small length scales, the criti- 
cal and fast-flow regimes typical of the RM (or domain 
wall) depinning, and at large length-scales, the critical 
and fast-flow regimes typical of the RP (or charge-density 
wave) depinning. From the equilibrium behavior of these 
kind of systems we have also inferred a richer dynamical 
roughness diagram including the low temperature creep 
regime which additionally includes roughness sectors cor- 
responding to the equilibrium geometry of the RP and 
RM classes. 

Our results are relevant for understanding the geome- 
try at depinning of periodic arrays of elastically coupled 
thin manifolds in a disordered medium such as driven 
particle chains or vortex-line planar arrays since these 
periodic systems display localized periodic correlation 
peaks. In particular our results are relevant for prop- 
erly controlling the effect of transverse periodic bound- 
ary conditions in large-scale simulations of constant-force 
driven disordered interfaces. From the roughness dia- 
grams of Figs. [2] and [10] we see indeed that the aspect 
ratio relation must be carefully chosen when taking the 
thermodynamic limit, depending whether one wants to 
study the large-scale behavior of a pure RM or a RP 
system. 

We have also argued that there is a geometrical equiv- 
alence between the d-dimensional periodic elastic system 
moving in d-dimensions and the d-dimensional elastic in- 
terface moving in a d-|-l dimensional periodic medium, 
although the mapping between these two systems is not 
exact. In this respect we note that the d = \ case we 
have studied numerically is the most stringent case, since 
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it goes beyond the usual small slope approximation used 
to develop the density in periodic components, Eq. HHl 
Indeed, since the roughness exponent for chains is larger 
than one, the average difference between the displace- 
ments of neighboring particles grows with the system size, 
violating the small slope approximation for large systems. 
Despite this fact, the results still remain valid even for 
the one dimensional case. We thus conclude that this 
equivalence is rather robust and should hold for higher 
dimensional cases. 
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Appendix: Elastic String vs Elastic Chains 

In this Appendix we show how the disorder Hamil- 
tonian for the elastic string in periodic disorder and for 
the periodic chain in one-dimensional disorder can be ob- 
tained. 



1. Two-dimensional periodic disorder 

The disorder contribution to the full Hamiltonian of 
the system is 



Hdis = J dx dzV{x, z)p{x, z). 



(A.l) 



The line density p{x, z) gives the position of the interface. 

In order to consider the periodicity of the system in 
the x-direction, the disorder potential can be written as 
a sum over periodic images as 



(A.2) 



where V{x, z) is the disorder potential defined in the in- 
terval X e [— Ajf/2, A//2]. In terms of a traditional un- 
correlated Gaussian disorder Vo(x, z), one can directly 
define 



V{x, z) = Qnix, z)Vo{x, z), 



using that 



Qm{x,z) = 



1 if x e [-M/2,M/2], 
otherwise. 



(A.3) 



(A.4) 



The disorder term thus becomes 

iJdis — J dx dzV{x, z)S{x ~ u{z)) 

= / dx dz 8j\f (x, z)Vq{x, z) 6{x — pM ~ u{z)). 

p 

(A.5) 

In order to obtain a disorder averaged Hamiltonian we 
use the replica-trick^^SiZ^ 

Hdis^-^Y. J dxdx'Vix)Vix')paix)pbix'). (A.6) 

ab 

Then, using the integral representation of the S function 

p{x, z) = S{x ~ u{z)) = -L /" dA e*^(^-"(^» , (A.7) 
2n J 

the replicated Hamiltonian reads^ 
ffdis = - ^ X! X! dxidzi dx2 dz2 

a.b Pi,P2 
X Om{xi,Zi)Qm{x2,Z2) 



X Vo{xi,Zi)Vo{x2,Z2) 

X S{xi ~ piM - Ua{zi))5{x2 -P2M - Ub{z2)). 

(A.g; 

Now, using for the disorder potential correlator that 



Vo{xi, zi)Voix2, Z2) = AoS{xi - x2)S{zi - Z2), (A. 9) 
and since 6^/ = 0, 1 implies 8^^ — Qm, one has 

^^dis = ~ -^X^X! f dxdzQmix^z) 

a,b pi,P2 

X S{x - PiM - Ua{z))d{x - P2M ~ Ub{z)). 

(A.IO) 

Now, we perform the sum over the localization function, 
resulting in 



J2 Six ~ PiM ~ u,{z)) = YiY. e^^^(--«"(^)\ (A.ll) 
where Ki = 27r??i/Af, and we used that 

E^"^*'-§E^('^-^0. (A.12) 

j K 

Thus, the disorder Hamiltonian can now be written as 
A 



2T]VP 



j dxdzQM{x,z) 



a,b Ki,K: 



X giKi{x-Ua(z))^iK2{x-Ub{z}) 



(A.13) 
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Now, using that 



dxe'^'^^-'^^^'' = M6 



(A.14) 



-A//2 



performing the integral over x and summing over K2 , one 
arrives at 



a,b K 



(A.15) 



2. One-dimensional periodic chain 



The positions of the particles are given by Xj = + 
Uj = ja + Uj, where x'j = ja is the nominal position in 
the unperturbed lattice, uj is the displacement and a is 
the average lattice space. In order to treat the model one 
should go from the Uj variables to a continuum formula- 
tion. Through this relabelling process, the decomposition 
of the density in this one dimensional problem leads to a 
density field 



where the last expression is valid for dxU ^ 1, if 
27rn/a, and the continuum field 



u{x) 



27r/a 



dq 
2^' 



iqx \ ^ ^if/a: 



3 



(A.17) 



is valid for x ^ a. This relabelling process is carefully 
described in Ref. [g^, especially Appendix A, and we refer 
the interested reader to this work. 

Considering the decomposition of the density, the 
Hamiltonian part corresponding to the uncorrelated 
Gaussian disorder Vq{x) is 



iJdis = / dxVQ{x)p{x) 



po / dxVo{x) 



(A.18) 



= Pa 



1 — dxu{x) + 



^ ^iK(x-u(x)) 



,(A.16) When replicating the Hamiltonian one has 



Hdis — — 



2T 



ab 



1 - dxu'^ix) ~ dxu\x) + ^^Kix-u^x)) ^ J2 e'K'{x-u\x)) ^ dxu''{x)dxu\x) 



- dxu^x) J2 e*^'(^-"'(-» - dxu\x) e^^(--""(-)) + ^ e--(^+^' 



)^-i{Ku'^{x)+K'u''(x)) 



(A.19) 



In this last expression one should drop constant shift terms and rapidly oscillating terms. Then, setting K = —K' in 
this last expression one has 



2T 



Yjdx 



dxu''{x)dxu\x) + Y 



-iK{u''(x)-u''(x)) 



(A.20) 
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